Nonextensive effects on the phase structure of QHD 
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We investigate nonextensive effects on phase transition in nuclear matter in the context Walecka 
many-body field theory. Changes in the general behaviour are observed when the results calcu- 
lated for the nuclear matter at vanishing baryon density are compared to those obtained through 
the standard Fermi-Dirac distribution. It is observed a dependence between the nonextensive pa- 
rameter q and the coupling constants C% of the phase transition. A numerical relation for this 
thermodynamical dependence is also proposed. 
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I. INTRODUCTION 

Non-extensive statistical mechanics (NESM) is a gen- 
eralization of thermodynamics and statistical mechan- 
ics aiming to overcome a number of physical systems 
that possess exotic properties, such as broken ergodic- 
ity, strong correlation between elements, multi fractality 
of phase-space and long-range interactions. The NESM 
framework proposed by Tsallis [l| is based on the non- 
additive g-entropy 

where fc is a positive constant, W is the number of mi- 
croscopic states, and pi is a normalized probability dis- 
tribution. 

In fact, additivity for two probabilistically independent 
subsystems A and B is generalized by the following pseu- 
doadditivity: 

S,{A,B) _ S,{A) S,{B) . S,{A)S,{B) 
fc fc ^ fc '^^ fc2 

For subsystems that have special probability correlations, 
extensivity is not valid for Boltzmann-Gibbs (BG) en- 
tropy, but may occur for Sq with a particular value of 
the index g ^ 1, called the g-entropic parameter. In the 
limit g — > 1, BG entropy is additive, i.e, 

Sbg{A,B) ^ Sbg{A) + Sbg{B) , 

and Eq. ([1]) reduces to the usual Boltzmann and Gibbs 
formula Si — k^^^pilnpi, which is the fundamental 
expression describing systems that usually do not present 
the exotic properties described above. 
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Several consequences of this generalized framework 
have been investigated in the literature (see also [| 
for a regularly updated bibliography) . This includes sys- 
tems of interest in high energy physics, e.g., the problem 
of solar neutrino 0, j^, the charm quark dynamics in 
a thermal quark-gluon plasma for the case of collisional 
equilibration Q, interpretations for central Au-Au colli- 
sions at RHIC energies in a Relativistic Diffusion Model 
(RDM) among others 

Recently, an interesting connection between quan- 
tum statistics and Tsallis NESM has been presented in 
Ref. [l6i]. The main result of this study is that, for g > 1, 
a g-generalized quantum distributions for fermions and 
bosons are given by 

nq{^l, T) = ^ ^^ I 1 ' (2) 

eq(/3(e - /i)) ± 1 

where reads 

r [l + (g-l)2:]^ if x>0 
~eq{x) ^ (3) 

[ [1 + (1 - g)a;]T^ if x<0. 

and X = P(e — fj,). In the g ^ 1 limit, the standard Fermi- 
Dirac distribution, n(fi,T), is recovered. As physically 
expected, as T — > 0, nq{p,T) —>■ n{fi,T). This amounts 
to saying that for studies of the interior of neutron stars 
(where, in nuclear scale, T ~ 0) we do not expect any 
nonextensive signature. On the other hand, in heavy ions 
collision experiments or in the interior of protoneutron 
stars, with typical stellar temperatures of several tens 
of MeV (1 MeV= 1.1065 x lO^" K), nonextensive effects 
may appear. 

In this paper, by following the arguments of Ref. [ll[, 
we discuss non-extensive effects on nuclear matter at zero 
baryon density. To this end we adopt the relativistic 
phenomenological approach developed by Walecka [l0|, 
the so-called quantum hadrodynamics (QHD). By taking 
into account the results on the phase transion in nuclear 
matter discussed in Ref. [l^, we show that there exist 
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FIG. 1: The solution M* /M of Eq. (O for nuclear matter (7 = 4) at vanishing baryon density as a function of temperature for 
the same value of C% and two different values of the parameter q. The right panel shows the same results but with a stretched 
temperature region around the point of phase transition. 



a thcrmodynamical dependence between the statistical 
correlations characterized by non-extensive parameter q 
and the coupling constant C|. A numerical relation in 
the q X C| plane resulting in regions of difFerents thcr- 
modynamical behaviors is also proposed. 



II. NONEXTENSIVITY ON QHD 

In the mean-field approach (for details see reference 
[To|). we have the scalar density 

es^^ I ^^Mv,T) + n,{~v,T)]dPk , (4) 



(27r)3 J E*{ky 



where E*{k) = V^^ + M^, /? = l/fcsT, M* is the ef- 
fective mass 



M* = M - g^cj = M~ ^ps 



(5) 



The parameter v = fi — g^^o = A* ^ [duj I'rriujY Qb is the 
effective chemical potential, and 7jv is the multiplicity 
factor (77V = 2 for pure neutron matter and 7Ar = 4 for 
nuclear matter). 

Since the function eq {x) is a deformed exponential, the 
mathematical convergence of the integrals in Eqs. ([l])-® 
must be verified. In Ref. [ll| it was pointed out that, to 
simultaneously satisfy the convergence of all integrals in 
Eqs.Q-®, we have found that 1 < g < 5/4, consistent 
with the limits of q in the range q e (0, 2] of Ref. [13] ■ 
Additionally, we also use for the coupling constants the 
values of reference namely [l9j. 



= 11.798 fm^ and ( 

m„ I \ m. 



8.653 fm^ , (9) 



which are fixed to give the bind energy -Ebind 
MeV and kp = 1.42 fm"\ 



-15.75 



and, instead of the usual Fermi-Dirac distribution, we 
have used for g > 1 the g-generalized quantum distribu- 
tion for fermions, 71,(1/, T), given by Eq. 

The baryon number density, the energy density and 
pressure are given, respectively, by 



Qb = 



7n 



(z/,T)-n,(-z.,T)]d3fc, (6) 



^ 9I 2,^9^ 
2 m^, 2 TO"^ 



7n 
(27r)= 



E*{k)[ng{iy,T) + ng{-iy,T)]d^k , (7) 



1 „2 lf'.(M-M*)2 + 



2ml^^ 2m2 
1 7N f 



3(27r)3 J E*{k) 



[nq{y, T) + nq{-v, T)]d^k , (8) 



III. THE PHASE STRUCTURE AT ZERO 
BARYON DENSITY 

In this Section, by considering the same arguments of 
Ref. ^15 ]. we explore the non-additive phase structure of 
the effective Lagrangian at vanishing chemical potential 
and baryon density {qb = v = Q). To this end, we first 
consider the range of values for the coupling constant C| 
given by 



(10) 



in the coupling-constant plane shown in Fig. (1) of 
Ref. [1^. At vanishing chemical potential, the terms with 
the baryon density do not appear in Eqs. ([7]) and In 
what follows, we will take as an example C| — 365. 

Theis et al. [lB| showed that the order of phase tran- 
sition is strongly dependent on the actual value of the 
coupling constant C|. However, in Sec. V of Ref. [ll| . 
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FIG. 2: The total energy density et and the scalar-field energy density to- divided by the q = 1 Stefan-Boltzmann limit 
(given by Eq. (26) of Ref. [ll|) as a function of temperature, at zero baryon density of nuclear matter (7 = 4). The same value 
of C| is considered for two different values of the parameter q. In the right panel, we show the same results in the stretched 
temperature region near the transition point. 




FIG. 3: The total pressure divided by the corresponding q — 1 Stefan-Bolztmann limit (Eq. (26) of Ref. 11]) as function of 
temperature. The right panel shows the first order phase transition point for q (solid curve) in the stretched temperature 
range. For q = 1.0218 (dotted curve) the phase transition is of second order. 



it is shown that we can avoid Maxwell construction by 
the variation of the parameter q. Since the order of 
transition depends on C| and the natural question 
is whether there exists some relation between C| and q 
at zero baryon density. 

In Fig. (1), the sudden drop in M* around T ~ 185 
MeV determines the abrupt rise of the energy density 
and pressure. We note that for q = 1 (Fermi-Dirac 
statistics) the self-consistency equation ([5]) has three so- 
lutions around T ~ 185 MeV imposing a sudden rise in 
the energy density and a peak in the specific heat. This 
behavior is shown in Figs. (2)-(4), where the temper- 
ature dependence of total energy density, pressure and 
specific heat divided by the corresponding high Stephan- 
Boltzmann temperature limit (for q = 1), given by 
Eq. (26) of Ref. [ll|, are displayed. The behavior of 
the curves for q = 1 characterizes a first order phase 
transition with the pressure curve crossing itself twice at 
T ~ 185 MeV. We observe that in this region the energy 
density and presure are also triple valued and that the 
specific heat is negative. The value of q = 1.0218 was 
obtained by requiring that energy density and pressure 
to be single valued characterizing a second order phase 
transition with non- negative specific heat. This allows to 



obtain a relation between C| and q explained below [20| . 

A. C| — q relation 

In order to obtain such a relation, in what follows we 
investigate the thermodynamical behavior of the nuclear 
matter for several values of the coupling constant C| in 
the q X C| plane. For each value of C5, the corresponding 
parameter q, for which the transition is of second order, 
is determined. Thus, a curve in the q x C| plane is ob- 
tained as shown in Fig. (5a). Below this curve the phase 
transitions are of the first order and above it the ther- 
modynamical behavior is smooth. Let us now show how 
the calculation is made via specific heat. 

Differently from the treatment discussed in Ref. p^ . 
the mathematical structure of the self-consistency equa- 
tion in our approach is not simple, so that the calculation 
must be done numerically. We observe that the specific 
heat calculated from Eq. ([7]) is hnear in dM* /dT. Thus, 
whenever there is a sudden fall in M*, we see a peak in 
the specific heat. By writing 

_ de _ de dM* 
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FIG. 4: The specific lieat of nuclear matter (7 = 4) divided by the corresponding q = 1 Stefan-Boltzmann limit (Eq. (26) of 
Ref. [ll|]) as function of temperature for the same value of C| and two different values of the parameter q. In the right panel, 
the same results in the stretched region around the phase transition point. For q — 1 and q — 1.0218 the phase transitions are, 
respectively, of the first and second order. 
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FIG. 5: Panel (a): Values of the critical parameter q for which the phase transitions are of second order as function of the 
coupling constant C|, for nuclear matter (7 = 4) at zero baryon density. Panel (b): The temperature corresponding to Panel 
(a) as function of qc- Panel (c): The same as in Panel (b) but for the effective mass. 



we can see from Eq. ([S]) that 

dT 1 + 2Cm* ^n,dfc ^ ^ 

where 

Cm- - {ga/m„fjN/n^ = Cl/M^. 

The singularities of dM* /dT lie in the curve determined 
by the vanishing of the denominator. Using Eq. ([5|), this 
condition becomes 

f°° dk 

M-2Cm'M*'J^ ,,^__ = 0. (13) 

Note that, for (7 = 1, we fully recover Eq. (18) of Ref. p^ . 

The number of intersections of the solutions obtained 
from Eq. ([T3|) and the self-consistency equation given by 
Eq. ^ determines how decoupling happens. We have a 
first or a second order phase transition, respectively, for 
two or one intersections. If there is no intersections, the 
thermodynamical behavior is continuous. The numerical 
results shown in Fig. (5a) can be summarized as follows: 



• For q lying below the q x C| curve the phase tran- 
sitions are of the first order. 

• For q lying on the q x C'g curve the phase transitions 
are of the second order. The corresponding values 
of temperatures and effective masses are shown in 
Figs. (5b) and (5c). 

• For q lying above the q x C| curve the decoupling 
is continuous. 

Therefore, the order of phase transitions depends not 
only on the actual values of C| but also on the values 
of q. 

IV. FINAL REMARKS 

We have investigated some effects of the Tsallis NESM 
on the mean field theory of Walecka (QHD) [10]. Instead 
of the standard Fermi-Dirac nucleon and antinucleon dis- 
tribution functions, we have used the g-quantum distri- 
bution obtained in Ref. [l^. Although there is no long 
range interactions in nuclear matter (Coulombian inter- 
action is not included in our analysis), it is worth empha- 
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sizing that the NESM effects appear due to the strong 
statistical correlations between elements of nuclear sys- 
tem. 

We have examined the phase structure of nuclear mat- 
ter at high temperature and at zero baryon density. The 
effective Lagrangian of QHD theory is considered for the 
same set of values of the coupling constants in the cou- 
pling constant-plane of Ref. Given that, at vanishing 



baryon density, the term with the baryon density do not 
appear in the equation of state, a relation between the 
coupling constant C| and the parameter q is obtained. 
As discussed, such a relation determines, in the q x C| 
plane, regions of different thermodynamical behaviors. 
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